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We consider the time evolution of simple quantum systems under the influence of random fluc- 
tuations of the control parameters. We show that when the parameters fluctuate sufficiently fast, 
there is a cancellation effect of the noise. We propose that such an effect could be experimentally 
observed by performing a simple experiment with trapped ions. As a byproduct of our analysis, 
we provide an explanation of the robustness against random perturbations of adiabatic population 
transfer techniques in atom optics. 
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I. INTRODUCTION 



In the last ten years the interest around the control 
and manipulation of quantum systems has grown very 
fast. The possibility to encode and process informa- 
tion has lead to innovative proposals. The major results 
have been achieved in Quantum Information Processing 
(QIP), including both theoretical and experimental ones. 
Quantum cryptography and information transfer pro- 
tocols |2J have enhanced our understanding of informa- 
tion processing and this in the near future will presum- 
ably lead to a significant technical advancement. Quan- 
tum Computation (QC) is still in a initial stage: even if 
the quantum computers seem to be able to solve quickly 
some problems which are intractable with classical com- 
puters |3j, more quantum algorithms are required to ex- 
tend its applicability. 

Unfortunately, the quantum system are very delicate 
and they are subject to two different kinds of errors. On 
one hand, there is the loss of information due to deco- 
herence the unavoidable interaction of quantum systems 
with their environments. This problem has been exten- 
sively studied over the past few years and proposals to 
overcome it have been put forward (and a few have been 
experimentally tested). These proposals include error 
avoiding error correcting strategies |5j and decou- 
pling techniques @. The other source of errors is the 
imprecise control of the parameters which perform the 
quantum operation (e.g. the laser or the magnetic field). 
How to handle such errors is an open problem, though 
some progress has been made in the framework of the so 
called geometrical quantum computation 0, H, Ej ■ The 
goal of this paper is to approach the second problem in 
very simple and idealized situations. 

A simple way to model the parameter noise is to con- 
sider a quantum system subject to a stochastic fluctu- 
ating field with zero mean. Such a model has been re- 
cently considered in order to study the effect of the noise 
on holonomic quantum gates: in Ref. [lfj it has been 
shown that there is a cancellation effect for a fast fluctu- 
ating stochastic field and shown that such a cancellation 
is due to the geometrical dependence of the holonomic 
operator. Recently, the general validity of such cancel- 



lation effect has been clarified: in Ref. [Tj is has be 
shown that for sufficient fast fluctuating stochastic field 
with zero mean the effects of the noise are wiped out. 

In this paper we shall study simple quantum systems 
subjected to stochastic noise and discuss some applica- 
tions. In Section ^ we consider random perturbations 
which are diagonal in the logical basis and propose an ex- 
periment to test the cancellation effect of the the noise — a 
simple modification of the experiment done by Kielpinski 
et al. 14]. Then we consider a more general noise and 
compute the fidelity. By elementary perturbation theory, 
we show that the effects of the noise are wiped out. In 
Section IIIII we discuss how noise cancellation could be 
relevant for adiabatic population transfer experiments in 
atoms optics (in this case, the cancellation effect avoids 
the break of the adiabatic approximation, and allows for 
desired transformation, even in presence of fast fluctuat- 
ing noise). 

For sake of simplicity, all the simulations and analytical 
calculation are done with Gaussian noise distributions 
but presumably analogous results can be achieved with 
a generic stochastic noise with zero mean |ll). 



II. TWO-LEVEL SYSTEMS 

Consider a two-level system evolving according to 
Hamiltonian 



H{t)=H + 8H!{t). 



(1) 



Suppose that Hq = B z o~ z (with o~ z being a Pauli matrix) 
and that SHj(t) is of the form 



N 



(2) 



where S A> are random 2x2 matrices and Sj (t) are "box 
functions" with time step r, i.e., they are functions equal 
to 1 in the time interval (Jt, (j + l)r) and zero otherwise. 
Let T be some "final" time at which we wish to consider 
the system, and let r = T/N; r can be regarded the 
correlation time of the noise. Then the evolution operator 



2 



from time zero to time T, generated by H (f) with SHj(t) 
given by @, can be written as 



where 



U(T) = U N U N - 1 .--U r -U 2 U 1 



Uj = cxp [-it (B z a z + 5A j )] 



A. Diagonal Noise 



(3) 
(4) 



First, we shall focus on the very simple case corre- 
sponding to the choice 5 A 3 — 5B 3 z a z , where the SB 3 ., 
j = 1,...,N are independent Gaussian distributed ran- 
dom variables with mean zero. In this case the noise 
is diagonal in the logical basis {| f), | J.)} in which a z 
is diagonal. Then J2J is trivially computed (due to the 
commutativity of the 8 A 3 ) , 



U(T) = U°(T)5U(T) 



(5) 



where U°(T) = exp(—iTHo) is the evolution generated 
by Hq, and 



5U{T) 








(G) 



A standard performance estimator is the fidelity |15| . 
which, in our case, is given by 



(^)(TM(T))\ 



(7) 



where ip^ Q \T) is the "ideal" final state evolved according 
to Hq (that is, when the noise is turned off) and ip(T) is 
actual final state (that is, when the noise is turned on). 
The two final states are generated, of course, by the same 
(generic ) initial state \tp(0)) = a\ T) + P\ !)• From © 
and JfjJ) one easily obtains 



T 2 (x, a, (3) = M 4 + \/3\ 4 + 2\a\ 2 \(3\ 2 cos x , (8) 



where 



1 



JV 



2T— V SB 3 
N ^ z 



(9) 



Since the random variables SB 3 Z in © are taken to be 
independent and identically distributed according to a 
Gaussian with zero mean and variance SB 2 , the proba- 
bility distribution of \ is 



Hx) 



N 



2tt2TSB 



_g 2(2T)' 2 6B' 2 



(10) 



By averaging JSJ with respect to IjlOl) . we arrive at the 
mean square fidelity: 



d X V( X )^T(x,a,P) 
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FIG. 1: The square fidelity from noisy process described by 
Q. The solid line is the theoretical mean square fidelity, eq. 
l(TTjl . for T = 100 (in arbitrary units) and SB = 0.1. 



As it can be easily seen, for N — > oo and r — > (while 
keeping Nt — T of order 1), the mean square fidelity 
T"^{a,P) approaches 1 (since \a\ 2 + \(3\ 2 = 1), indepen- 
dently of the initial state. This corresponds to a short 
correlation time of the noise, that is, to a fast random 
fluctuating field. 

Now we'd like to provide some perspective on such 
limit behavior. Of course, the limit r — > is only an ide- 
alization for t small but finite. The above results are ob- 
tained for constant SB 2 , but the Heisenberg uncertainty 
relation imposes strong constraint on the energy fluctua- 
tion happening in such a short time interval, SB 2 oc 1/r. 
Note that for SB 2 oc N eq. {TDJ and {TTJ describe a 
system interacting with a white noise environment that 
is the standard way to model the decoherence effect. So, 
for very small r the variance SB 2 should not be any more 
considered constant and our approximation of constant 
variance breaks down. However, rather than microscopic 
(environmental) noise, we are interested in modeling the 
macroscopic (parametric) noise due to imprecision in the 
control field, and for such our approximation should ap- 
ply. In this case it is interesting to ask whether the con- 
dition of small t could be physically relevant. In general, 
in a quantum evolution the final time T is fixed; the cor- 
relation time r can be hardly controlled (i.e., stabilizing 
the control field) and we are not in the condition to have 
a cancellation effect. In Section IIIBI we present an ex- 
perimental proposal to test the presence of this effect; 
the experiment realizability lies in the control of the cor- 
relation time t of the simulated noise. Moreover, there 
are situations in which this effect can be experimentally 
relevant: when we have a further degree of freedom and 
can change the evolution time T. In these cases, fixed 
t with the above properties, we can prolong T in order 
to have T ^> r (that is N ^> 1) and exploit the cancel- 
lation effect. We give an example of such situation in a 
adiabatic evolution in Section ITTT1 

We observe that another interesting feature of eq. I|ll|) 
is that the cancellation of the noise does not depend on 
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FIG. 2: The decoherence effect during the evolution for two 
different value of N = T jr. Theoretical curves (dashed line) 
and those obtained by numerical simulations (solid line) are 
showed. For the theoretical curve (eq. I13H the decoherence 
time is tdeco = 2t 2 5b 2 ' 



the strength SB: it is always possible to find a suitable 
./V in order to obtain cancellation of the noise — invert 
equation I jllj l and express TV in terms of SB for given Tt 
and evolution time T. 

We recall that dependence of the fidelity on a specific 
choice of the initial state is usually eliminated by aver- 
aging !Fj,(a,/3) over all the possible initial states with 
respect to the uniform distribution on the unit sphere in 
the Hilbert space of the system in our case the (pro- 
jective) sphere \a\ 2 + |/3| 2 = 1. Performing this operation 
yields to 

^=<^(a J /3)>=i(2 + e- 32 ^ si ) (12) 

The simple model we have been considering here is 
often used as a toy model for phenomenological decoher- 
ence 0|. The relationship between noise cancella- 
tion and decoherence is easily seen by considering the 
time evolution at the discrete times t = kr, k < N. By 
proceeding as before, the average square fidelity at time 
t = kr results 

— o 1 / 2TtSB 2 , 

^ 2 = -(2 + e »-) (13) 

with which it is quite natural to associate a 'decoherence 
time' t deco = 2T N gB± (see Figure |2 Qj, [l3| ) . 

B. Experimental proposal 

In Ref. 0| Kielpinski et al. used the same idea to 
simulate the decoherence effect due to interaction of the 
quantum system with the environment degree of free- 
dom. The authors used trapped ions and study the 
coherence of superposition of quantum state subject to 
simulated noise. The logical states were the hyperfine 
states of a trapped Beryllium ion \F = 2, m,p = — 2) and 



\F = l.mj? = —1) sublevels of the ground state 2 Si/2- 
The environment noise was simulated shining the ions 
with a off-resonant laser with random varying amplitude 
for the electromagnetic field E 3 and random intensity 
(proportional to (E 3 ) 2 ). The laser electromagnetic field 
produce a AC Stark effect on the ions and let one state 
to acquire a random phase respect to the other. This 
effect is quadratic (quadratic Stark effect) in the electro- 
magnetic field E 3 . In fact, the two hyperfine states have 
the same angular momentum (they have both L = 0) 
and the difference is in the spin part of the wave func- 
tion. The splitting of the energy level is linear (linear 
Stark effect) for state with different angular momentum 
since only these states have non-vanishing matrix ele- 
ment (L = i\E ■ r\L = k) (with i ^= k). Then in the 
above example we have corrections to the energy level 
only quadratic and not linear in E 3 . Because of the 
quadratic energy shift, we have a random phase differ- 
ence proportional to (E 3 ) 2 for every r interval. Even 
if this effect is sufficient to produce decoherence effect 
(as found by the authors), we cannot presumably see the 
cancellation effect discussed above in a transparent way 
since our new stochastic variable \ x ^2(E 3 ) 2 nas 110 
zero mean. 

A small modification of this experiment should allow 
us to see sharply this cancellation effect. To have an evo- 
lution described by Hamiltonian Q it is sufficient to use 
states with different angular momentum in order to pro- 
duce a linear Stark effect. At every time interval t the 
perturbation of the laser produce a shift of the energy 
levels proportional to the random intensity of the laser 
E 3 ; this produces an evolution where the phase difference 
between the states is given by a known dynamical part 
plus a random phase exp(2iaE 3 r) (where a is a propor- 
tional constant). In this case, the new stochastic variable 
X oc J2j E 3 has still zero mean and we expect to obtain 
results shown in the previous section: fixed the evolution 
time T we should see an increase of T 2 as the corre- 
lation time decreases (see Figure^. Moreover, if the 
environment decoherence does not depend on the simu- 
lated noise, we should be able to see effect analogous to 
the ones in Figure In particular, subtracted the effect 
of the environment, the decoherence time should increase 

&S tdeco 0^ ~ ■ 

C. Off-Diagonal Noise 

We now consider the case of off-diagonal noise, that 
is, the matrices SA 3 in (J2J) are of the form SA 3 — 
5B 3 X <7 X + 5B 3 y (jy, with SB 3 ,, SB 3 y , j = l,..,,N indepen- 
dent Gaussian distributed random variables with mean 
zero, and a x , <j y being the usual Pauli matrices. Then 
the one-step evolution operators Uj in © are 

Uj = e - igJ - 3T = cos(B 3 t) - in 3 ■ asm(B 3 r) , 

where B 3 is the modulus of the vector B 3 = 
(SB X ,SB 3 \ B z ) and n 3 is the associated unit vector. By 
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FIG. 3: Cancellation effect for system subject to off-diagonal 
noise discussed in section 111 CI For the stochastic process 
SB = 0.01. When the number of random extraction increases 
the fidelity approaches 1. 



FIG. 4: Figure shows the population evolution of the three 
states during the adiabatic transfer process for two different 
correlation time noise (ji — T/Ni). For the first curve iVi = 
10 4 and for the second one N 2 = 10 5 for T = 100 (in arbitrary 
units) . 



a perturbation expansion in = 8B k /B z , k — x,y, we 
obtain 

Uj = cos(B z t) — ia z sm(B z r) 

ARJ SB 3 

Note that the zero order terms are nothing but U°(t). 
We shall denote by SUj the first order terms, i.e., 



5BI 



SBl . 



5Uj = i(a x "—2-+ a y — ^)sin(£ z r) 



Under the assumption that e/. = 8B J k /B z <C 1 the 
operator U (T) can be easily computed by first order per- 
turbation theory: by keeping in © only the terms that 
are first order in we obtain 



U(T) = U°(T)-J2 6p3 



(14) 



By taking into account that sin(i? z T) w B z t and 
exp(iB z r) w 1 for r = T/N < 1, we may finally evalu- 
ate the modulus of the scalar product and compute the 
fidelity. We obtain 



1 - 2TRe 



N N ' y 



It is important 
N — > oo and r 



(15) 

to note that, also in this case, for 
(while keeping Nt = T of order 
1), the fidelity approaches 1. This is so because 8B k are 
independent random variables with mean zero [Isj . 

So, also in this case there is a noise cancellation ef- 
fect. This effect is confirmed by the numerical simula- 
tions shown in Figure 



III. ADIABATIC EVOLUTION 



where U°(T) = cxp(-iTH a ) and 

5Pl = {U {r)) j ~ X 8U 3 (U a {T)) N - 3 

To compute the fidelity 10 for the (generic) initial 
state |^(0)) — a\ |) + (3\ I) we need to calculate the 
scalar product (ip^(T)\t/j(T)) = (U°(T)ip(0)\U(T)ip(0)). 
According to (|14|l . this is given by 

l-<^(0)|(C/ (T))t(^^)|^(0)) 



The matrix elements of 8P J in the logical basis are: 
CLjSP^mm = and 



CE Sp3 > 



Sm(B z T) y-> pl B^T(N-2j+l] 



B 7 



(8Bi±i5Bi). 



The adiabatic population transfer is an important tech- 
nique used in atoms optics to achieve population trans- 
fer between quantum states of atoms and molecules . 
We first create coherence between the initial and final 
state (population trapping) and then produce an adia- 
batic evolution to transfer the population to the final 
state. This scheme has seen a great success and has 
been used in many different areas: chemical reaction pflj . 
laser-induced cooling atoms optics , cavity quan- 
tum electrodynamics |23l |24} . The wide range of applica- 
tion is due to many advantages of this scheme: it is easy 
to implement in different system, it has an high rate of 
population transfer and it is robust respect to variations 
of field parameters [25|. 

Consider two states |1) and |2) coupled to an excited 
state |e) by two lasers (i.e., a A system). The states 
|1) and 1 2) can be degenerate or quasi-degenerate but 
it is important that we can address separately both of 
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FIG. 5: Figure shows the fidelity and the population of exited 
state |e) as function of T/t = TV. 



them. The Hamiltonian in the rotating frame with reso- 
nant laser frequencies is 



H 



-(ni(i)|l)<e|+fi 2 (*)|2)<e|)+/i.& (16) 



where £li(t) and ^(i) are the time-dependent Rabi fre- 
quencies and depend on the parameters of the lasers (am- 
plitude and phase) . The diagonalization of l|16|) gives two 
eigenstates \B±) = -^(±fi|e) + Q^l) + Q 2 \2) (called 
bright states) respectively with eigenvalues ±f2(t) = 

±\jY!n=i l^il 2 ; an d an eigenstate \D) = 1/0(0 2 |1) — 
Oil 2)) (called dark state) with zero eigenvalue. In the 
adiabatic evolution (i.e. when the JVs change slo wly 
and QT 3> 1) it follows from the adiabatic theorem |2y| 
that if the system starts at time t = in an eigenstate 
of H(0) (dark or bright state) during all the evolution 
it will remain in the eigenstate of H(t) with the same 
eigenvalue. 

Now we provide a simple example of the foregoing. 
Suppose that f2i(i) and ^(t) are such that fl is time- 
independent and Sli (0) = and ^2(0) = ft. Moreover, 
suppose that the initial state is |1). Then slowly turn 
on the first Rabi frequency and turn off the second one. 
The system will always be in the dark state \D) and, 
at the end of the evolution (i.e., when fii(T) = fi and 
0,2 (T) — 0), we will be in |2) state and have achieved 
population transfer from |1) to |2). 

Consider now the case in which we have not a com- 
plete control of the laser field but the Rabi frequencies 
can fluctuate fii — > f2i + Sfli (where the 5Qi are inde- 
pendent and Gaussian distributed with zero mean and 
variance a 2 ). This effect can produce errors in popu- 
lation transfer scheme for two reason: in general noisy 
perturbations may yield to significantly different output 
state and fast fluctuations could break the adiabatic ap- 
proximation leading to transition to undesired (bright) 
states. 

Numerical simulations show that, again, in the fast 
fluctuation regime the noise effects cancel out. In Fig- 



ure 01 we present the population evolution of |1), |2) and 
|e) states subject to noisy evolution during the |1) — > |2) 
coherent adiabatic transfer |28j. More precisely, we start 
from |1) and, during the evolution, the |2) is populated; 
at the end only the |2) is present. The |e) is never popu- 
lated because of the high value of f2T parameter. These 
simulations are done using QT = 1000. Since we are in 
the adiabatic regime we are sure that the errors present 
are those induced by the perturbations. In Figure 0| we 
show the population evolution when the system is subject 
to noise with different correlation time r (and different 
N = T/t). For TV = 10 4 the transfer operation is not pre- 
cise (i.e. the |2) state is not completely populated) and 
the \e) state is populated; this is a sign of the breaking of 
the adiabatic approximation due to noise. For TV = 10 5 
the evolution is much more similar to the ideal one: |2) 
is completely populated and |e) never appears during the 
evolution. 

The above results can be explained calculating the am- 
plitude transition from the dark to the bright state in 
presence of perturbations. The standard rule (see, e.g. 
27]) to calculate the probability amplitude of a transi- 
tion from the n state to one of the k state (k ^ n) at time 
t is a fe = E'nlloo^^expiif^knt^dt' where £' 
indicates that the term n — k is omitted, u>kn — Ek — E n 
and Vkn is the matrix element associate to the transition 
n — > k. In our case, the initial state is the dark state 
and the final states are the bright states. The eigenvalues 
±fl are constant and then LOkn = fi- The perturbation 
Hamiltonian for time jr < t < (j + 1)t in the |e), |1), 

|2) basis is SV^ = J2i=i SCl^\i){e\ + h.c. and in the new 
(dark-bright states) basis the relevant matrix element are 



V DB, 



-V 3 
v DB- 



3 



S(jT, (j + l)r) 



dV k , 



(17) 

we must take 
= S(t - jr) - 



To calculate the matrix element —^j- 
into account that dS(jr,(j + l)r)/dt 
S(t — (j + 1)t). Let us focus our attention only on 
the first term in l|17|l . by differentiating it we obtain 
l/(v / 2fi) J2j SQi(dn 1 /dt)S(jT, (j + l)r) - fii<5fi J 2 ((5(t - 
jr) — S(t — (j + 1)t). The terms d£li/dt representing the 
adiabatic driven evolution are very small and can be ne- 
glected. Inserting this result in the expression for and 
performing the integration we have 



fl fe = ~^Yl m 2i a n(jT)Ol(jT)e 



iQjT 



- a n ((j + l)r)Qx((j + l)T)e m ^ T + 

+ terms with (Oi «-> O2) (18) 

If a n and f2j change slowly i.e., a n ((j + l)r) w cL n (jr) 
and fij((j + 1)t) w fii(j'r). The exponential terms can 
be simplified to obtain exp(iilr(j + 1/2)) sin(fir/2)/(2i) 



6 



and eq. JTHJ for r = T/N < 1 gives 



T 



a k 



e l ° JT « r 



(jr) 



n r- ^2 



Iv" 

(19) 

The sums of converge to the mean of <5f2j that is 
to zero; the other factors are bounded (0 < \a n \ < 1 and 
< | fij | < Q) and by consideration similar to those at 
the end of section III CI (and pij ) we can conclude that 
afc — > for N — > oo. In our case, |n) = |_D), |fc) = |£>±) 
and |a_B ± | 2 — > as iV - >oo : the transition from dark to 
bright states is suppressed in the fast fluctuating regime 
and the evolution happen in the dark space. 

To have a more detailed picture in Figure |3] we show 
also the trend of the fidelity (upper curve) and the rela- 
tive average population of the |e) state (lower curve) as 
functions of T/t. The trends of the two curves are corre- 
lated, which is suggesting that the main source of error 
in the operation is the population of excited state due 
to loss of the adiabatic approximation. As expected, be- 
cause of the cancellation effect, for great T/t the fidelity 
approaches 1 even in presence of strong noise and the |e) 
state is not populated. 

These results not only can explain why the adiabatic 
population transfer scheme is robust against field fluctua- 
tion but can give information for the experimental set-up. 
In fact the adiabatic evolution is, in general, arbitrarily 
long; once the experimental parameters are fixed (as the 
laser with its proper noise correlation time r) we can pro- 
long the evolution time in order to increase T/t and let 
the noise average out. As shown before, for every noise 
strength a and correlation time r we can find and evolu- 



tion time T in order to obtain the desired fidelity: that 
is, to achieve the population transfer with arbitrary small 
error. 



IV. CONCLUSION 

We studied the effect of stochastic noise on several 
quantum systems. The noise is described by a Gaussian 
stochastic process with zero mean superposed to the ideal 
quantum evolution. For all of the systems we found that 
for fast fluctuating noise a cancellation effect appears: 
the noise fluctuations average out leading the system to 
a state near to the ideal one. 

We showed by analytical and numerical calculation 
how this effect can appear in a two level system and 
propose an experiment to verify the presence of this can- 
cellation regime. The experiment is based on the one 
performed in Ref. 0] and we think that, with a mod- 
ification of the experimental set-up, it could be easily 
performed. 

We applied the same model to another important tech- 
nique in atoms physics : the adiabatic population trans- 
fer. We explained how this effect leads to the robustness 
of the adiabatic process against the perturbation noise. 
This can be important for the experiments using the adi- 
abatic population transfer. 



Acknowledgments 



P.S. wishes to thank D. Kielpinski for useful comments. 



[1] C.H. Bennett and G. Brassard Proceedings of IEEE In- 
ternational Conference on Computers, Systems and Sig- 
nal Processing 175-179, IEEE, New York, 1984. C. H. 
Bennett Phys. Rev. Lett. 68(21), 3121, 1992. 

[2] C.H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. 
Peres, and W. K. Wootters Phys. Rev. Lett. 70, 1895, 
1993. 

[3] P.W. Shor Proceeding of 35th Annual symposium on 
foundation of Computer Science. (IEEE Computer So- 
ciety Press, Los Alamitos, CA, 1994). L.K. Grover Proc. 
28 th Annual ACM Symposium on the Theory of Compu- 
tation, 212, ACM Press, New York, 1996. 

[4] P. Zanardi and M. Rasetti, Phys. Rev. Lett. 79, 3306 
(1997). 

[5] P.W. Shor Phys. Rev.A 52, 2493 (1995); A.M. Steane 
Phys. Rev. Lett. 77, 793 (1996); E. Knill, R. Laflamme, 
Phys. Rev.A 55, 900 (1997) and references therein. 

[6] L. Viola and S. Lloyd, Phys. Rev. A. 58, 2733 (1998); 
L. Viola, E. Knill, and S. Lloyd Phys.Rev.Lett. 82, 2417 
(1999); D. Vitali, P. Tombesi, Phys. Rev. A 65, 012305 
(2002); P. Zanardi, Phys. Lett. A 258 77 (1999). 

[7] J. A. Jones et al. , Nature 403, 869 (2000). G. Falci et al. , 
Nature 407, 355 (2000). P. Zanardi and M. Rasetti, Phys. 



Lett. A 264, 94 (1999). G. Falci et al. , Nature 407, 355 

(2000) . R.G. Unanyan, B.W. Shore and K. Bergmann, 
Phys. Rev. A 59, 2910 (1999). L.-M. Duan,J. I. Cirac 
and P. Zoller, Science 292, 1695 (2001). L. Faoro, J. Siew- 
ert and R. Fazio, Phys. Rev. Lett. 90, 028301 (2003). I. 
Fuentes-Guridi et al. Phys. Rev. A 66, 022102 (2002). A. 
Recati et al. Phys. Rev. A 66, 032309 (2002). P. Solinas 
et al. Phys. Rev. B 67, 121307 (2003). P. Solinas et al. 
Phys. Rev. A 67, 062315 (2003) 

[8] WangXiang-Bin, M. Keiji Phys. Rev. Lett. 87, 097901 

(2001) ; WangXiang-Bin, M. Keiji Phys. Rev. Lett. 88, 
179901(E) (2002). X.-Q. Li et al. Phys. Rev. A 66, 042320 

(2002) . S. L. Zhu, Z.D. Wang, Phys. Rev. Lett. 89, 
097902 (2002). PSolinas et al. Phys. Rev. A 67, 052309 
(2003) 

[9] A. Blais and A.-M. S. Tremblay Phys. Rev. A 67, 012308 

(2003) ; A. Nazir, T. P. Spiller, and W. J. Munro, Phys. 
Rev. A 65, 042303 (2002); G. De Chiara, G.M. Palma, 
Phys. Rev. Lett. 91, 090404 (2003); A. Carollo et al, 
Phys. Rev. Lett. 90, 160402 (2003); A. Carollo et al 
quant-ph/0306178 V.I. Kuvshinov, A.V. Kuzmin, Phys . 
Lett. A, 316, 391 (2003); F. Gaitan, quant- ph/0312008 



7 



[10] P. Solinas, P. Zanardi, N. Zanghi quant-ph/0312109 

[11] P. Facchi, S. Montangero, R. Fazio, S. Pascazio 
quant-ph/0407098 

[12] Our sampling set of the Bloch sphere (i.e. our inizial state 
space) is given by {±ei} i=Xt y )Z , (±e x ± e v )/s/2, (±e x ± 
e z )/y/2, (±e y ± e z )/\/2. Here d denotes the normalized 
vector of the i-th direction. 

[13] We make simulations with 1000 different realizations and 
average the resulting square fidelity. A part from a noisier 
in the curve in the plots, the results are weakly affected 
by this parameters so that with 100 realizations we obtain 
the almost the same results. 

[14] D. Kielpinski et al. Science 291, 1013, (2001). 

[15] M.A. Nielsen and l.L. Chuang Quantum computation and 
information, Cambridge University Press (2000). 

[16] D.M Meekhof et al, Phys. Rev. Lett. 76 1796, 1996; 
S. Schneider and G.J. Milburn Phys. Rev. A 57, 3748, 
(1998); M. Murao and P.L. Knight Phys. Rev. A 58, 663 
(1998). 

[17] M.O. Scully and M.S. Zubairy Quantum Optics, Cam- 
bridge University Press (1997). 

[18] This is an elementary instance of the law of large num- 
bers: Suppose to have a sum written in the form 
Ylf=o f( T j)~iv where Xj are independent random variable 
with zero mean and variance a 2 , r = T/N and < j < T. 



Var(YN) (the variance of Yn = jg> 52j=o f( T j) x i) i s 

W E^Lo(/( r J')) 2 - Thus > if / is bounded in the inter- 
val [0, T], Var(Vjv) — > when TV — > oo and therefore 
Yn — > since the Xj mean is zero. Since in our case 
f(rj) = exp(—i2B z rj) the second terms in eq. I15H goes 
to zero and T — > 1 for N — > oo for every initial state (i.e. 
independently from a and 0). 

[19] K. Bergmann, H. Theuer, and B.W. Shore Rev. Mod. 
Phys. 70, 1003 (1998). 

[20] P. Dittmann et al. J. Chem. Phys. 97, 9472-9475 (1992). 

[21] S. Kulin et al. Phys. Rev. Lett. 78, 4815-4188 (1997). 

[22] M. Weitz, B. C. Young, and S. Chu Phys. Rev A 50, 
2438-2444, (1994). 

[23] A.S. Parkins et al. Phys. Rev. Lett. 71, 3095-3102, (1993) 

[24] R. Walser, J.I. Cirac, and P. Zoller Phys. Rev. Lett. 77, 
2658-2661, (1996). 

[25] L.P. Yatsenko, G. Guerin, and H.R. Jauslin 
|quant-p h/0100706^. 

[26] A. Messiah Mecamque Quantique, Dunod Paris (1964). 

[27] L.I. Schiff Quantum Mechanics, McGrawn-Hill (1968). 

[28] We are interested in the |1) — > |2) transfer then we start 
form |1) and do not sample the initial state space. The 
fidelity is still obtained averaging the results of 1000 dif- 
ferent simulations. 



